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Abstract. A review and critique of the Isgur-Paton flux tube model of hadronic physics is presented. 
I> \ This entails a detailed comparison with recent lattice gauge theory results which exposes both the 

successes and shortcomings of the model. Applications to hybrid masses, meson and hybrid meson 
^ | ' decays, hadronic charge radii, the spin-orbit force, baryonic hybrids, and hybrid photocouplings are 

also discussed. Finally, I comment on the issue of adiabatic surface crossing which appears in both 
the flux tube model and lattice studies. 

c5o FLUX TUBES 

<N ■ 

The impending upgrade of CEBAF at Jefferson Lab and the advent of a new experi- 
! mental facility, Hall D, devoted in large part to the exploration of the gluonic excitations 

of QCD make it a propitious time to review the flux tube model. The flux tube model 
(FTM) of Isgur and Paton is now nearly 20 years oldJH] and much has been learned in 
the meantime. This note shall review the foundations and predictions of the FTM and 
compare it, where possible, to lattice gauge or other theoretical tests. Extensions of the 
original ideas and new applications are also treated. 

The flux tube model is an attempt to form a tractable model of low energy QCD which 
is based on longstanding ideas of string-like gluonic degrees of freedom. The central 
idea is that gluons should rearrange themselves into flux tubes which, in the heavy 
quark limit, adjust their configuration instantaneously in response to quark motion. Thus 
quarks are constrained to move on adiabatic gluonic energy surfaces. The lowest such 
surface is the familiar 'Coulomb+linear' potential of the constituent quark model and 
lattice gauge theory. States based on this surface thus form the 'conventional' mesons 
of the quark model. Higher adiabatic surfaces represent gluonic excitations and mesons 
built on these surfaces correspond to 'hybrids' [2]. Thus the flux tube model is a simple 
and intuitive way to extend the nonrelativistic constituent quark model to include gluonic 
degrees of freedom. 

The model was originally motivated through a series of truncations on the Euclidean 
time strong coupling QCD Hamiltonian with a lattice regulator (Hamiltonian lattice 
gauge theory). The degrees of freedom are quarks fields on lattice sites and gluonic 
'link variables' Ui = exp(—iagA^(x)) where £ represents the link (jc, fx). In the strong 
coupling limit the Hamiltonian is given by 
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where g is the strong coupling, a is the lattice spacing, and n is a lattice site. The velocity 
variables Ui have been replaced by electric field operators E(. Gauge invariant pure glue 
states are formed by closed (possibly multiply connected) loops of link operators. The 
commutation relation [E a ,U] = T a U then implies that the energy of these states is simply 
the sum of the quadratic colour charges of each link: 

£loop = f- 52 
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where ^ 2 = 4/3 for a field in the 3 or 3 representations, 10/3 for 6 or 6, etc. The 
presence of quarks permits gauge invariant states with open flux strings which terminate 
on quark colour sources or sinks. Perturbations to these states are provided by subleading 
quark hopping and magnetic terms. The former allow flux tube breaking via quark pair 
production or quark motion. The latter can change link colour representations, cause link 
hopping, or change loop topology. 

Isgur and Paton simplify the dynamics by (i) assuming an adiabatic separation 
of quark and gluon degrees of freedom (ii) neglecting 'topological mixing' such as 
loop breaking or loop Euler number changing transitions (iii) working in the nonrel- 
ativistic limit. The model is meant to be applied to the 'intermediate regime' where 
l/a ~ \fb. They then model link variable dynamics in terms of spinless colourless par- 
ticles ('beads') of mass ba where b is the string tension in the static quark potential. 
Finally these particles are assumed to interact via a linear potential and perform small 
oscillations about their resting positions. The result is a simple discrete string model for 
glue described by the Hamiltonian: 
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where y n is the transverse displacement of the nth of N string masses , p n is its momen- 
tum, bo is a bare string tension, and R = (N + \ )a is the separation between the static 
quarks. This Hamiltonian may be diagonalised in the usual way yielding 

Hftm = boR + {^R - ~ - ^ + • • •) + £ a nX (4) 

where at, creates a phonon in the nth mode with polarization A . Notice that the string 
tension has been renormalized by the first term in brackets. The last term in brackets 
is the Liischer term of string phenomenology Q3|] . The mode energies are given by 
(O n = 2/as'm[nn/2(N + 1)]. Thus (0\ — > n/R as N — > °° is the splitting between the 
ground state Coulomb+linear potential and the first gluonic excitation surface at long 
range. The energy of a given phonon state is approximately 
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N = m ( n m++n m -) (6) 
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where n m ± is the number of left (right) handed phonons in the mth mode. 



HYBRIDS 

Hybrid mesons are constructed by specifying the gluonic states via phonon operators 
and combining these with quark operators with a Wigner rotation matrix: 
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The projection of the total angular momentum on the qq axis is denoted by A = 
T.m( n m+ — n m -). The parity and charge parity of these states are given by 

P\LM L ;SM s ;A,{n m +,n m -}) = (-) L+A+l \LM L ;SM s ; -A, {n m _,n ra+ }>, (8) 

C\LM L ;SM s ;A,{n m+ ,n m -}) = (-) L+s+A+N \LM L ;SM s ;-A,{n m -,n m+ }). (9) 
States of good parity are thus formed as 

\LM L ; SM S ; £; A, {n m+ , n m _ }) = (\LM L ; SM S ; A, {n m+ , n m -}) + £ \LM L ;SM S ; -A, , n m +})) 

(10) 

Possible single phonon (m = 1) mesons are listed in Table I. Underlined quantum 
numbers represent quantum number exotic hybrids (mesons with quantum numbers not 
available to a qq state). 

TABLE 1. Some Single Phonon 
Mesons. 
g L S J PC 

+10 1++ 

+ 1 1 (2,1,0)+- 
+ 2 2~ 

+ 2 1 (3,2,1)-+ 



1 1" 

1 1 (2,1,0)-+ 

2 2++ 

2 1 (3,2,1)+" 



Isgur and Paton obtained hybrid meson masses by solving a model Hamiltonian of 
quark motion on the single-phonon excited surface: 

2/1 dr L 2jlr z 3r r 



The interaction term incorporates several important additional assumptions. Namely the 
njr phonon splitting is softened at short range. The parameter / which appears in the 
softening function was estimated to be roughly unity[l]. Furthermore, it was assumed 
that the attractive Coulomb (1/r) potential remains valid for hybrid mesons. Note that 
this is at odds with the expected 

Vo g e = +f (12) 

potential of perturbative one gluon exchange (the colour factor is that appropriate to 
gluon exchange between quarks in a colour octet state). This is an important assumption 
which will be examined in more detail below. 

Finally the quark angular momentum operator is now complicated by the presence of 
gluonic/string degrees of freedom. One may write 



Lqq =L — Ls^ —Ls ± (13) 

where L (L5) is the total (string) angular momentum. Note that L$ L mixes adiabatic 
surfaces. Using Lc. = Ar and neglecting surface mixing yields the centrifugal term of 
Eq.EU This additional assumption will also be examined in the next section. 

The hybrid masses obtained by solving Eq.[TT]are labelled Ejp in Table II. Isgur and 
Paton also estimated the effects of adiabatic surface mixing and used these as their final 
mass estimates (labelled E' IP ). The column labelled KW is explained in the next section. 
Finally, Table I implies that hybrids with quantum numbers 2 ±T ,1 ±T ,0 ±T ,1 ±± are all 
degenerate at this order in the FTM. 

TABLE 2. Hybrid Mass Predic- 
tions 



flavour 


Ejp 


F' 


Ekw\4\ 


1=1 


1.67 


1.9 


1.85 


1=0 


1.67 


1.9 


1.85 


ss 


1.91 


2.1 


2.07 


cc 


4.19 


4.3 


4.34 


bb 


10.79 


10.8 


10.85 



The next section will present a survey of possible tests of this 'zeroth order' Flux Tube 
Model. 



TESTING THE FLUX TUBE MODEL 

Small Oscillation and Adiabatic Approximations 

The small oscillation approximation may be tested by considering a model of trans- 
verse beads interacting via a linear potential. Numerically solving such a Hamiltonian||5|] 
reveals that the small oscillation approximation is accurate for long strings but overesti- 
mates gluonic energies by an increasing amount as the interquark distance shrinks. Typi- 
cal energy differences are order 100 MeV at 1 fm. Similarly, the adiabatic approximation 
can be tested by numerically solving the coupled quark-bead system. One finds||5|] that 



the adiabatic approximation underestimates true energies by roughly 100 MeV, with 
slow improvement as the quarks get very massive. It thus appears that these approxima- 
tion errors tend to cancel each other, leaving the IP predictions intact. 

Gluonic Surfaces 

Recent advances in computational speed and algorithms make it possible to test some 
of the assumptions of the FTM against the predictions of lattice gauge theory. Figure 
[l] shows the assumed ground state Coulomb+linear (solid line) and first excited state 
(dashed line) potentials of Isgur and Paton. These are compared to lattice Wilson loop 
computations of the same interactions (points)||g]. It is apparent that the IP potential 
overestimates the strength of the Coulomb potential and underestimates the string ten- 
sion. Of course both of these model parameters are obtained by fitting the meson spec- 
trum and therefore include effects which are not in the Wilson loop. 
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FIGURE 1. Comparison of the Potentials of Eq.^^with lattice gauge theory|6|. 



Should the Coulomb Potential Appear? 

More troubling is the first excited potential (bursts) which shows signs of saturating 
(or perhaps turning repulsive) at small distances, in disagreement with the assumption of 
Isgur and Paton. Indeed, simply omitting the attractive Coulomb portion of the IP hybrid 
potential dramatically improves agreement of the model with the lattice. 

This disagreement represents something of a conundrum because Isgur and Paton 
had a good physical reason to employ the colour singlet qq Coulomb interaction in 
their model: if one assumes the repulsive short range interaction of Eq. [121 it becomes 
energetically favourable for the system to emit a gluon once the interquark separation 
becomes small enough. This gluon combines with the 'valence' gluon of the hybrid 
to form a scalar glueball, thereby changing the quark colour configuration to that of a 



singlet, and the Coulomb potential to the attractive form of Eq. [TT] The point at which 
this should happen is indicated by the arrow in Fig. El It is clear that the lattice sees no 
such behaviour. This may be simply because the minimal relative momentum permitted 
on the lattices employed in the study were too coarse to permit the decay. 

A physical reason for the suppression of this coupling is also possible. For example, 
if one considers the hybrid to be dominated by a Fock space component consisting of 
a constituent quark, antiquark, and gluon, then the postulated transition occurs through 
gluon emission from either the valence gluon or a valence quark. In the former case 
the coupling to a glueball is zero due to the colour overlap while in the latter case the 
coupling is suppressed by the large (infinite) quark mass. Thus one expects the coupling 
between the surfaces to be very small, which implies that the surface mixing will not be 
seen unless lattices with exceptionally large temporal extents are employed. 

Which short distance behaviour should the model use? It seems clear that surface 
mixing which changes Fock sectors should not be used in a potential model. Rather, 
such mixing should be incorporated by explicitly including the appropriate transition 
operator in the formalism. Thus, for example, flux tube breaking is expected to eliminate 
the linear potential for distances beyond roughly one fermi when light quarks are present. 
However, such a truncated linear potential should not be used to construct mesons, rather 
one should employ the linear potential to form a basis of bound states and then include 
mixing to four quark states in an appropriate form. This approach avoids the unpleasant 
situation of having a mesonic ionization energy. Similarly, the hybrid potential should 
not include the mixing term since this mixes hybrid with hybrid+glueball Fock states. 




FIGURE 2. 
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Adiabatic Surface Crossing. 



Higher Surfaces 

The FTM predicts an entire tower of hybrid surfaces and it is interesting to examine 
their form in light of lattice data. Juge, Kuti, and Morningstar have carried out a detailed 
analysis of the relationship of the hybrid surfaces of Fig. 1 to the string excitations of 



Eq. [5J7] . They have found that surface excitations only have n/r splittings for very 
large source separation (roughly 4 fermi or greater). This is shown in Fig. |3] where the 
dashed lines are the predicted Nn/r energy differences of the FTM. This is something 
of a surprise since one expects a phonon-like excitation spectrum on general grounds. 
It appears that QCD strings are complex objects at intermediate distance scales. Finally, 
the figure shows a cross over region at about 1 fermi where the surfaces move from a 
perturbative behaviour (characterized by the 'gluelump' spectrum) to a more string-like 
behaviour. 
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FIGURE 3. Hybrid Surface Energy Differences. 2ro is roughly 1 fermi|7]. 



Hybrid Masses Revisited 



In light of these issues it is appropriate to revisit the original hybrid mass estimates 
of Table I. A number of variants of the IP calculation were made in Ref. [4]. Before 
presenting some of these, note that there is additional ambiguity in the quark angular 
momentum term of Eq. HH namely it is also possible to set L q q = L — J g where J g is the 
total gluonic angular momentum. Squaring then yields 

L 2 qq = L(L+l)-2A 2 + (J 2 g ). (14) 

It has been found that setting J g = 2 yields good agreement with lattice data lEfljl Cthis 
was also found in a constituent quark model 112 ill ). Notice that this is not numerically the 
same as the centrifugal term in Eq.fTTI 



Employing the hybrid lattice potential of Fig. 1 and Eq.UH yields the hybrid masses 
labelled Ej^w m Table U. It is seen that these can differ by up to 200 MeV from the 
IP predictions of the second column. The similarity to the third column is a fluke since 
adiabatic surface mixing is not included in this computation. Merlin and Paton have 
noted that the majority of adiabatic mixing effects may be absorbed into the static 
potential by including the moment of inertia of the string in the centrifugal term (see 
the discussion below): 



2jur 2 2^2 + 1^3' 

and with a more important effect which modifies the strength of the njr splitting to 
be larger as r becomes larger than m q /b. Thus the final revisited prediction for a light 
hybrid is roughly 2.1 GeV. 



THE IKP DECAY MODEL 

Shortly after its introduction, the flux tube model of meson structure was extended by 
Isgur, Kokoski, and Paton to provide a description of meson[8] and hybridll] decays. 
The transition operator was envisioned as arising due to the quark hopping term of the 
lattice QCD Hamiltonian. The lowest terms in the expansion of this operator are 

Hhop = £ vla-iiy n + a£ v4« 1 M v HYn- (16) 

If one assumes a smooth string then the first term dominates as the lattice spacing gets 
small and one has a 3 5i strong decay operator. Alternatively, if the string is rough then 
the first term averages to zero upon summing over all local string orientations and the 
second term dominates, yielding a 3 Pq strong decay operator. The authors of Ref. 
assume the second scenario since it is supported by experiment lllOll . 

Flux tube degrees of freedom were incorporated by assuming factorization: 

({. . .}bd;{. . .}bd\0\{. . .}bU f ) w (bd;bd\ 3 P \bW) ■ ({. . .}; {. . .}|{. ..}) (17) 

The first matrix element on the right hand side is a typical 3 i 5 o mesonic decay overlap. 
The second represents the overlap of the gluonic/flux tube degrees of freedom. Assuming 
that the quark pair creation occurs at a transverse distance y± from the interquark axis 
of the parent meson yields the results 

({0...};{0...}|{0...}>~e-M (ig) 

for meson decay and 

({0...};{0...}|{l...}}~ J± e-/^ (19) 

for hybrid decay. The factor / is a computable constant of order unity. The extra factor 
of v_j_ in the hybrid decay vertex forces the decay to pairs of identical S- wave mesons to 



be zero. This is the origin of the famous 'S+D' selection rule in this model (it occurs in 
other models for different reasons). 

Unfortunately, it will be some time before lattice gauge theory has progressed to 
the point where models such as these can be thoroughly tested. However, preliminary 
computations reveal that a substantial closed flavour decay mode (such as bbg — > Xb 7 !) 
may exist[ll]. In the meantime we hope that experiment will provide some clues. An 
alternative decay model is discussed in the next section. 



The hadronic decay model of the proceeding section is its most well known extension. 
However, the model has been extended in several other directions as well; some of these 
are described here. 



Several years ago Isgur pointed out that the energy carried by the flux tube will change 
several features of the naive quark model[12] (see also jllfll V For example, zero point 
oscillation of the flux tube about the interquark axis will induce transverse fluctuations 
in the quark positions, something which is not present when the flux tube is treated as 
potential. The additional fluctuations have the effect of increasing the charge radius of a 
heavy-light meson (qQ): 



where the second term in the bracket is the new contribution. He estimated this to give 
rise to a 50% increase in charge radii of light quark hadrons. 

This observation has subsequently been expanded upon by Close and Dudek who 
observe that radiative decays of hybrid mesons may proceed because the recoil of the 
radiating quark affects the string degrees of freedom giving a nonzero overlap of the flux 
tube wavefunction with the ground state flux tube wavef unctions of ordinary mesons. 
Preliminary computations with this mechanism have appeared ! 1 311 . 

A similar scheme involving the emission of pointlike pions may be used to compute 
hybrid decays to final states such as ^p lfLUl . The most striking result here is that this 
decay mechanism evades the 'S+D' suppression discussed above. 



EXTENSIONS 



Charge Radii and Hybrid Decays 




(20) 



Adiabatic Surface Mixing 



Merlin and Paton examined the effects of adiabatic surface mixing on the leading 
order FTM by considering the complete quark-bead system ab initio a 15II . Although the 
effects can be quite complicated, with mixing to all surfaces possible, they found that 



the majority of the effects can be absorbed in a redefinition of the hybrid potential by 
including the rigid body moment of inertia of the string in the centrifugal term and by 
modifying the strength of the ic/r term (see Eq.fT3l. 

An explicit computation revealed found mass shifts of order -100 MeV for conven- 
tional S-wave light quark mesons and +200 MeV for light quark hybrids. The resulting 
mass splittings are quoted in column three of Table I and were used by Isgur and Paton 
to form their final estimates of the lowest lying hybrid masses. 



Spin Orbit Forces I 

Merlin and Paton also examined spin orbit forces in the context of the FTMlfl6ll. The 
idea was to map the operators of the leading spin orbit term in the heavy quark expansion 
of QCD, namely Vso = g/2mo- B, onto FTM degrees of freedom (phonons). Merlin and 
Paton did this by identifying the magnetic field with the lattice operator 

T a B a {x)~^(u P {x)-Ul{x)). (21) 

Since the plaquette operator moves flux links in a fixed topological sector, it is natural 
to identify the magnetic field with the bead kinetic energy. Doing so then allows one to 
write Vso in terms of phonons. 

Explicit computations revealed that spin orbit splittings due to Vso are small and that 
the majority of the splittings arise from Thomas precession, Vjh = \(jq x r q ) ■ a. This is 
modelled by including the effects of phonons on the quark coordinate (see the discussion 
of the charge radii above). The mass splittings for light hybrids are listed in Table III. 
One sees that the lowest member of the octet of light hybrids is predicted to be the 2 + 
while the heaviest is the + ~. This appears to be in conflict with lattice gauge theory 
which finds that the lightest hybrids are 1 h . 

TABLE 3. Spin Orbit Hybrid Mass Splittingsd . 

J PC | 2+~ 2~+ 1-+ Q-+ 1+- 0+- 1++ 1 — 

(5M(MeV) I -140 -20 20 40 140 280 



Spin Orbit Forces II 

Spin-dependent forces in the FTM were also taken up by the authors of Ref. ifnll in 
an attempt to resolve a conundrum in the spin-orbit sector of the quark interaction. The 
issue is that many models of hadrons prefer a vector Dirac structure of confinement, 
rather than the phenomenologically accepted Dirac scalar interaction. This may be 
studied by using the heavy quark Foldy-Wouthyusen version of the QCD Hamiltonian 
in Coulomb gauge. The resulting &{\jm) and &{\/m 2 ) operators, H\ and H2, depend 
on the chromoelectric and magnetic fields and hence are nonperturbative. In a similar 
approach to that of Merlin and Paton, the authors of Ref. 11711 chose to study these 
operators by mapping the chromo fields to phonons. In this case the mapping was based 



on the idea that the electric field counts string length and therefore should be mapped as 

(x = no) = -^p (y^ (n + 1) — y a ^ (nj) where a is a colour index and A is a polarization 
index. Imposing the canonical commutation relations gives the magnetic field as a 
derivative with respect to transverse displacement. Finally both expressions can be 
mapped to phonon degrees of freedom with the standard Fourier transform. Notice that 
this approach has the flux tube beads carrying colour charge. 

The form of the effective spin-dependent interactions were then evaluated by inserting 
these field operators into H\ and H2. It was found that an effective scalar spin orbit 
interaction could indeed arise in the flux tube picture of hadronic structure even though 
the static confining interaction was vector. This was the case because nonperturbative 
mixing of mesons with hybrids contribute to spin-dependent interactions and can change 
the naive expectations based on the nonrelativistic reduction of a simple interaction 
kernel. 



Vector Decay Model 

The same chromofield-phonon mapping which was used in the study of the spin orbit 
interaction 11 711 was used to examine hybrid decays in Refs. Ill 8ll . In this case the operator 
of interest is simply V = —g J i/AOC • Ayr. The resulting decay vertex is given by 

2 1 

H in*= lJ ^L I ^COs( W 0^(^)ff-«A(^fifi)(<*-<!l)4^ r Qfl)' (22) 



where the e(r) are polarization vectors orthogonal to f . The integral is defined along the 
QQ axis only. This model also has an 'S+D' decay rule, however in contrast to the IKP 
decay model this is due to a node along the interquark axis (the cos(^^)) which causes 
the selection rule, rather than a node perpendicular to the interquark axis. 



The phenomenology of this model has been presented in the second of Ref. Ill 811 . One 
finds that it is similar to model of meson decays in that D wave decay modes tend 
to be suppressed. Recent comparison with experiment have proven surprisingly accurate 
and lend support to hybrid interpretations for nonexotic 2 h (2003) and l ++ (2096) 
states |ET 



Hybrid Baryons 



Capstick and Page have made a detailed study of baryon flux tube dynamics II22I1 . This 
is a technically challenging problem due to the multitude of vibrational and rotational 
modes which are available to a Y-shaped string system. However, they have found 
that the problem simplifies considerably because the string junction decouples to good 
accuracy from the rest of the bead motion. Thus a hybrid baryon can be approximated 
by three quarks coupled via linear potentials to a massive 'junction bead' . The dynamics 
of this system are completely specified by the FTM and variational calculations indicate 

that the lowest lying hybrids are J p = j + and | + states at approximately 1870 MeV. 



Happily, lattice investigations of the static baryon interaction have begun l23ll . The 
chief point of interest is whether the expected flux tubes form into a 'Y' shape or a 'A' 
shape. This may be addressed by carefully examining the baryonic energy in a variety 
of quark configurations. Current results are mixed, with some groups claiming support 
for the two-body hypothesis [24] and some for the three-body hypothesis [25]. Finally, 
a strong operator dependence in the flux tube profiles has been observed ll26ll . which 
clearly needs to be settled before definitive conclusions can be reached. 

CONCLUSIONS 

The flux tube model is now nearly 20 years old. In this time it has been applied to 
an increasing array of problems and extended in several directions, chiefly by taking 
seriously the idea that dynamical string-like gluonic degrees of freedom are important 
in the low lying mesons. A number of additional extensions are: 

(i) glueballs (glue loops). A preliminary study has been made in the original paperjl]], 
however much remains to be done here. Comparison with lattice gauge theory will 
provide a crucial test. 

(ii) FT effects in baryons. The charge radii ideas of Isgur, Close, and Dudek have 
immediate impact on baryons and should be studied. 

(iii) adiabatic surfaces. It is worthwhile to attempt to leverage new precise lattice data 
on the gluelump spectrum and the hybrid adiabatic surfaces to improve the FTM in 
detail. 

(iv) the FTM may allow one to improve the semiclassical fragmentation formalismLlJ. 

(v) long range spin- spin and spin-orbit forces should be re-examined in an attempt to 
pin down this difficult aspect of nonperturbative QCD. 

In summary, the FTM provides a compelling picture of strong QCD dynamics; how- 
ever, it is a picture only ! We have seen that the FTM correctly describes the level or- 
derings and, perhaps, splittings of gluonic adiabatic energy surfaces at large distances 
(perhaps as large as 4 fm). The model fails to describe the spectrum at small interquark 
distances (although, of course, it can be amended). Furthermore, the original IP model 
of hybrids is likely to be incorrect in many details, although its phenomenology may be 
surprisingly robust. We have seen that it is possible to extend the model in many dif- 
ferent ways. Of course these extensions rely on detailed aspects of the FTM which are 
untested at best. It appears likely, for example, that the spin orbit splitting of Ref. Illfll 
do not agree with lattice data. In the end, the utility of a tractable and appealing model 
should not be underestimated. 
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